Self-consistent superconducting proximity effect at the quantum spin Hall egde 
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We study self-consistently a microscopic interface between a quantum spin Hall insulator (QSHI) and a 
superconductor (SC), focusing on properties related to Majorana fermion creation. For an s-wave SC we show 
that odd-in-momentum, or p-wave, order parameters exist for all doping levels of the QSHI and that they can be 
related to different spinless Cooper pair amplitudes. Despite this, the induced superconducting gap in the QSHI 
always retains its s-wave character, validating the commonly used effective model for superconductivity in a 
topological insulator. For a d xy -wa\e SC, we show that a Majorana mode is only created at finite doping and 
that there is no excitation gap protecting this mode. 

PACS numbers: 74.45,+c, 71.10.Pm, 74.90.+n 



Topological insulators (TIs) are a new class of materials 
where the bulk is gapped but there exist conducting surface 
states which are robust against all time-reversal invariant per- 
turbations. With the experimental discovery of a two dimen- 
sional (2D) TI, or equivalently, a quantum spin Hall insulator 
(QSHI) III], in a HgTe/CdTe heterostructure J2[] and the subse- 
quent advent of 3D TIs ft^, these materials have lately gained 
a large amount of attention |0, IH . One of the areas attracting 
considerable interest is when the TI surface state is supercon- 
ducting. Fu and Kane [6] showed that the low-energy spec- 
trum in a 3D TI with an s-wave superconducting order param- 
eter at high doping levels resembles that of a spinless p x + ip y 
superconductor (SC) and thus supports Majorana fermions at 
vortices [7] or at the interface between superconducting and 
ferromagnetic regions 10, [H. The non-Abelian statistics of a 
Majorana fermion can provide fault tolerant topological quan- 
tum computation [9] and multiple proposals already exist on 
how to detect Majorana fermions in TIs lioll . However, de- 
spite the far reaching consequences of the p-wave character 
of the superconducting state, very little is known about its mi- 
croscopic origins. For example, the Majorana mode survives 
even when the doping level goes to zero 10, 0] although there 
has so far not been any reported signs of p x + ip y -character 
in this doping regime. Moreover, superconductivity is in gen- 
eral only added on a phenomenological level as a constant 
order parameter to the effective Hamiltonian of the surface 
TI state |0, [H [Tot] despite being generated by proximity to a 
SC. Very recently Stanescu et al. [11] studied proximity in- 
duced superconductivity in a 3D TI and found, using several 
approximations, a p-wave order parameter only at high dop- 
ing levels. They also predicted that this component will even- 
tually close the superconducting gap, which would, in fact, 
render the commonly used phenomenological model with a 
constant s-wave order parameter invalid in the high doping 
regime. A comprehensive picture of the symmetries of the 
induced superconducting state in a TI is thus clearly lacking 
and, as a consequence, the connection to the existence of Ma- 
jorana fermions is not well established. 

In this Rapid Communication we will treat a QSHI-SC in- 
terface fully self-consistently, in which we do not only calcu- 
late the superconducting proximity effect, or the leakage of 



Cooper pairs into the QSHI, as has been done before ill ill , but 
also the accompanied loss of Cooper pairs, or inverse proxim- 
ity effect, on the SC side of the interface. Even more impor- 
tantly, we will explicitly calculate all significant p-wave am- 
plitudes and show that there are, in fact, large contributions 
even at zero doping. Our results demonstrate that it is natural 
that the Majorana mode survives even at zero doping. More- 
over, we show that the superconducting gap never closes but 
always retains its s-wave character which provides a formal 
motivation for the commonly used phenomenological model 

IM llOfl ■ We have also self-consistently studied the interface 
between a QSHI and a ci^-wave SC, since a d^-wave state in 
a TI was recently shown to support Majorana fermions II 1 211 . 
We show, however, that in our microscopic model the Majo- 
rana mode only appears at finite doping and what there is no 
excitation gap protecting this mode, thus rendering it useless 
for quantum computation |9j LL3|] . 

To self-consistently study the interface between a SC and 
the edge of a QSHI we consider a minimal microscopic model 
defined as H — H\ + Hs + and schematically displayed in 
Fig. Q] For H\ we use the first QSHI prototype, the so-called 
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FIG. 1: (Color online) Microscopic details of the QSHI-SC interface 
for the (1,0) surface of the SC and with sites a; a (dark), bi a (light). 



Kane-Mele model 11411 . which is defined on a honeycomb lat- 
tice with spin-orbit coupling A: H\ = -tJ2(i.j), a a \a a ja + 

^/E l!Q aL a « + * A E<( iJ )) ^j a L s a/3«j> Here a is the 
fermion operator on the bipartite honeycomb lattice, (i, j) and 
denote nearest neighbors and next nearest neighbors 
respectively, a, (3 are the spin indices, and Vij = +1 (— 1) 
if the electron makes a left (right) turn to get to the second 
bond, t is the nearest neighbor hopping amplitude and we 
set t = 1 for simplicity. We also fix A = 0.3 which gives a 
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bulk band gap of 1 and thus allows the chemical potential //i 
to vary between and 1 . We further define the SC on a 2D 
square lattice at half-filling with an on-site Hubbard attraction 
U to generate a prototype conventional s-wave superconduc- 

tor: H S = -tE<^>,a b V>i,« " ^E, We treat 

iJs within mean-field theory using the self-consistency con- 
dition = — U (bi^bif), for the s-wave order parameter 
Au- We will use U = 1.2 which gives A u = 0.1, but we have 
also verified our results for Ajj = 1. Finally, the tunneling 
Hamiltonian is described by — — ?E(j j) a a la^',a+H.c., 
where a is on the edge site of the QSHI, II, and b is on 
the first SC site, SI. t is technically allowed to vary from 
to 1, although we expect an experimental interface to never 
reach t = 1. We align the TI so that its edge is along the 
zigzag direction, connect the (1,0) surface of the SC to the 
QSHI, and assume the same unit cell size a for both the QSHI 
and the SC to avoid any lattice mismatch. We also assume a 
smooth interface and Fourier transform in the direction along 
the QSHI-SC interface. We can solve the above model self- 
consistently by first diagonalizing the Hamiltonian for a guess 
profile of the order parameter A;j(i) throughout the whole 
structure. Ajj(i) can then be recalculated by using the mean- 
field self-consistency condition. By repeating this process un- 
til Ajj(i) does not change between two subsequent iterations, 
we achieve self-consistency for the superconducting state in 
the whole QSHI-SC structure. The superconducting proxim- 
ity effect is explicitly captured in the (on-site) Cooper pair 
amplitude Fjj(i) = (a^a^) in the QSHI and, equivalently, 
with a — > b in the SC, where obviously also Fjj = —Ajj /U. 

Self-consistent superconducting gap. — InFig.|2ja) we show 
the self-consistent profile of Fu across the QSHI-SC interface 
for a few representative values of pi and t. As seen, there is 
a depletion of Cooper pairs, or inverse proximity effect, on 
the SC side of the junction with the accompanied leakage of 
pairs, or proximity effect, into the QSHI. The spread in Fjj 
on the QSHI side is essentially the same as the spread of the 
QSHI edge state, whereas the depletion region is a few unit 
cells wide on the SC side. As can be expected, the proxim- 
ity effect is larger for larger t, whereas, perhaps more sur- 
prisingly, there is very little dependence on the QSHI doping 
level /ii for small to moderately large tunneling amplitudes. 
Only a self-consistent approach can capture the inverse prox- 
imity effect and in Fig. 0b) we plot the value of Ay on the 
first SC site, SI, as a function of p,j. As seen, this value is 
significantly lower than Ajj = E^bulk) (dashed line). We 
also in (b) plot the extracted superconducting gap in the QSHI 
edge state, E g (ll), The self-consistent values of E g are shown 
with a solid line whereas non-self-consistent results, using a 
step-function of Ajj = _E s (bulk) at the interface, are shown 
with crosses and, as seen, they are essentially identical. There- 
fore, self-consistency is not important for calculating the in- 
duced superconducting gap at the edge of a QSHI from an 
s-wave SC. This is in contrast to a d xy -wa.ve SC where self- 
consistency is crucial as we will discuss below. We also see 
that E g (ll), in general, increases with pi for all experimen- 
tally relevant t. The slight downturn at small pi is due to 




FIG. 2: (Color online) (a): Fu across the QSHI-SC interface for the 
combinations (pi,t) = (0,0.2) in (black x), (0,0.6) in (blue □), 
(0.75, 0.2) in (red o), and (0.75, 0.6) in (green o). Vertical dashed 
line marks the interface and the distance is measured in number of 
unit cells. Inset: Fu as function of E g at site II for < fj,i < 1, 
t < 0.6, and the bulk order parameter Au = 0.1 (black) and 1 
(red). Crosses mark data, lines are linear fits, (b): — Aj/(S1) and 
i5 9 (II) as function of pi for t = 0.2 (black) and 0.6 (red). Crosses 
mark E g (ll) reached in a non-self-consistent calculation. Dashed 
line indicates \Au\ = E g in the SC bulk. 

the fact that for a highly transparent interface the SC induces 
a slight electron doping in the QSHI edge state whereas the 
(positive) applied chemical potential pi creates an overall hole 
doping. We note especially that even for the highest pi we can 
achieve and still have a bulk gap, i.e. /i] < 1, E g {\\) is always 
finite. This is in contradiction to the prediction by Stanescu et 
al. ifTHl for a 3D TI. While p-wave components are induced 
at the interface as we show below, we attribute the non-zero 
Eg to the fact that the s-wave component Fjj on site II is al- 
ways finite for any fij, and thus E g is set by this value of Fu- 
In fact, in the inset in Fig. |2ja) we show on a linear relation 
between Fjj and E g on site II, where the linear relation is in- 
dependent of both pi and t. To describe processes along the 
QSHI edge, it is therefore possible to use a model where the 
effect of the SC is only taken into account through an effective 
Aj; = U e sFu for the edge state, where the effective pairing 
potential U e g is independent of both pi and t. Such an effec- 
tive model has been completely dominating in the literature 
SIEOj and our calculations provide a fully self-consistent 
evidence for the validity of this approach. 

Induced p-wave order parameters. — As established by Fu 
and Kane [6], the low-energy spectrum in a 3D TI with an 
s-wave superconducting order parameter resembles that of a 
spinless p x + ip^-wave SC in momentum space, although 
time-reversal symmetry is not broken. However, in the present 
case of a 2D QSHI with a ID edge, any odd-in-momentum 
component necessarily only has a simple p-wave symmetry. 
Even though this p-wave symmetry is intimately linked to the 
generation of Majorana fermions, no detailed study has yet 
appeared on its origins. Stanescu et al. Ill 111 very recently pre- 
dicted, after some rather restrictive approximations, that the 
p-wave component is directly proportional to (lit 2 , but no de- 
tailed calculations were preformed. Here we explicitly calcu- 
late both on-site and bond superconducting order parameters 
that are odd in momentum, or p-wave. Since the s-wave state 
in the SC is induced through a local pairing and the QSHI edge 
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state is extremely localized at the interface, all significant p- 
wave pairing amplitudes should be captured within these two, 
locally defined, quantities. 

We start by defining the momentum-resolved on-site pair- 
ing amplitude Fjjk(i) = (ai-ki a ik^)- The conventional 
on-site order parameter can then be written as Fjj(i) — 
j^- J2k=-Tr ^Uk{i)> where N k is the number of fc -points. Ob- 
viously, only the part of Fjj k that is even in fc contributes to 
Fjj. There can also, technically, be a term in Fjj k that is odd 
in fc, which we will call F^ k . For A = this odd component is 
zero, but on the edge of a QSHI it is finite. We can understand 
this odd component better by defining the spinless on-site op- 
erator Cfc = -^(cik^+aki) at the QSHI edge. Then the Cooper 

pair amplitude (c fc c_ fc ) = \{a ki a^ - a_/.j.a fet ) = Fg k , 
due to the explicit coupling of spin and momentum at the edge 
of a QSHI. Thus, F{j = J2l=o F Uk is nothing else than 
a spinless on-site p-wave pairing amplitude. In Fig. Qa) we 
show the spatial extent of Ffj for a prototypical interface. We 
see that Ffi peaks at site II but that it is also present on site 
SI, and that there is a sign change between the two sites. This 
sign change interestingly also encodes an odd character across 
the interface, and we speculate that this might be a ID man- 
ifestation of the p x + ipy-wave present on the surface of 3D 
TIs. Fig. [3jb) furthermore shows that Fg increases approxi- 
mately linearly with fij, with t determining the slope. Since 
Ffj changes sign with fjti, the finite charge transfer between 
the QSHI and the SC discussed above makes Fg slightly nega- 
tive for small positive (bulk) ji\ at, especially, large t. Without 
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FIG. 3: (Color online) (a): Fu (thick black), F v (thin black), F° JX 
(thin red), and Fj 2 = F% (thin dashed red) across the QSHI-SC 
interface for pci = t = 0.5. Vertical dashed line marks the inter- 
face and the distance is measured in number of unit cells. Inset: 
fc-dependence at site II for F uk (black) and Fj lk (red), (b): Fu 
(thick) and F u (thin) on site II and Fj (thin dashed) on the t-bond 
for t — 0.2 (black) and 0.6 (red) as a function of fix. 

this charge transfer Fy is zero at /ii = 0. This can also be seen 
from the inset in Fig.[3fa), where we plot the fc-dependence of 
the p-wave order parameters. F^ k peaks sharply at the k value 
where we find the QSHI edge Dirac cone in the Bogoliubov 
spectrum. Since this cone is at fc = ±ir when /ii = 0, time re- 
versal symmetry will then force Ffi to vanish J^lOj]]. We also 
note that for all small to moderately large tunneling ampli- 
tudes Fy is linear in t, which is also different from the earlier 
prediction of P behavior 111 ill . We also note that F[} < Fjj for 
all hi and t. This further supports our earlier conclusion that 



the induced p-wave state will never close the superconducting 
gap in the QSHI, even at high doping levels. 

We also study the Cooper pair amplitude on nearest neigh- 
bor bonds, Fjs(i) = ■^(aiiai + sf-aifa i+ si), where 5 labels 
the bonds as defined in Fig.Q] A small, but finite, Fjg exists at 
the QSHI-SC interface, but its value is largely independent on 
A. Thus this (even) order parameter is not a consequence of 
the topological nature of the QSHI but it is a property of the 
surface of the SC. However, by using the same momentum- 
resolved technique presented above for Fjj, we can also define 
an odd-in-momentum component F° ]S , which is only non-zero 
for finite A. We find that especially Fj^Il), which resides on 
the t bond, is large. This could have been expected, since 
the QSHI edge state leaks substantially into the SC. Further- 
more, we find that F° ]S increases approximately linearly with 
t but that it is essentially independent on as seen in Fig. [3] 
Previous work JfjUTl]] has established a finite p-wave order pa- 
rameter at high doping levels, but our results explicitly show 
that there also exists a finite induced p-wave order parameter 
at [i\ = 0. This is despite time-reversal symmetry dictating 
that F° JSk = at k = ±vr for half-filling. The finite F° JS is 
still possible due to a rather broad fc-dependence, as seen in 
the inset in Fig. Oa). To further analyze Fj S we define the 
spinless bond operator d k = |[(a^ - b kj ) - {a ki - b H )}, 
where a is on the II site and b is on the SI site. Then we can 
construct the spinless Cooper pair amplitude residing on the t 
bond as (d k d- k ) = \ [i^(Il) + F$(S1)-V2F° 1 (II)] = F° } . 
Since F^ at sites II and SI have different signs, Fj^Il) con- 
tributes crucially to this spinless p-wave bond order parame- 
ter. In Fig. [3jb) wee see that Fj increases moderately with t, 
similar to F\j and Fjj, but, contrary to Fg, it is essentially in- 
dependent on /ii in the whole range < fij < 1. In summary, 
we have here explicitly shown that spinless p-wave order pa- 
rameters exists for all doping levels of the QSHI. At low to 
moderate doping levels the spinless bond order parameter Fj 
dominates, whereas for very high doping levels the spinless 
on-site Fy eventually becomes larger. This helps explain why 
the Majorana mode found in e.g. vortex cores exists for all 
doping levels of the TI. 

d X y-wave superconductor. — While the p-wave character of 
an induced s-wave superconducting state at a TI surface is 
of wide interest due to Majorana fermion creation, so are of 
course other systems with Majorana fermions. Recently Lin- 
der et al. IU2I1 discovered that Majorana fermions also exists 
at the surface of a d^-wave superconducting 3D TI. The sur- 
face/interface of a d^-wave SC contains zero energy states, 
or midgap Andreev bound states, since A(9) = —A(tt — 6), 
where 9 is the angle of incidence (see e.g. Ref. fl5l) . When 
these states exist on the surface of a TI, their double spin- 
degeneracy is broken and a Majorana fermion appears instead. 
This is contrast to the topologically trivial high-T c SCs where 
spin-degeneracy is never broken but gives rise to "double" 
Majorana fermions. From an experimental perspective, how- 
ever, there are two crucially important questions that needs to 
be addressed before realizing this Majorana system. (1): In a 
self-consistent model of a QSHI - d xy SC interface, will the 
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induced superconducting state in the QSHI be of the required 
d X y-wave character, and (2): how protected is the Majorana 
fermion, i.e. how large is the excitation gap? 

By rotating the SC in Fig. [T]to a (1,1) interface and, in- 
stead of U, using an effective pairing amplitude J on near- 
est neighbor bonds, we can model an interface between a 
QSHI and a d^-wave SC. We start the self-consistency loop 
with Ajs — —JFjs having a pure d^-wave symmetry, i.e. 
Aji = — Aj2, but we do not impose any symmetry con- 
straints on Ajs during the loop, as to not prevent the forma- 
tion of other, competing pairing symmetries at the interface. 
This procedure results in a solution that is fully d xy -wave in 
the bulk and persists out to the SI site with only a moder- 
ately large inverse proximity effect. This d^y-wave state also 
leaks weakly into the QSHI insulator where it very rapidly 
decays. However, there is a larger, but imaginary, compo- 
nent of Fji on the i-bond, accompanied by an imaginary and 
opposite signed Fj2 = Fj^. In aggregate this gives a local 
d xy + id x 2_ y 2-wave symmetry iflal at the QSHI edge. Thus 
an effective model with a purely d^y-wave superconducting 
order parameter at the QSHI edge is not valid for a prox- 
imity induced d^y-wave state. However, a small, but finite, 
dj;j,-wave component still exists at the interface, so zero en- 
ergy states might still exist at the interface. We thus turn to 
investigating the Bogoliubov band structure in Fig. |U The 




FIG. 4: (Color online) Bogoliubov band structure for the QSHI - d xy 
SC structure with [ij = t = 0.5 and J = 1. Unaffected by the 
interface are the Dirac cone on left hand side of the QSHI (red, close 
to k — — 7r) and Andreev bound states at the right hand side of the 
SC (green). Andreev bound states at the QSHI-SC interface (blue, 
crosses) (resolution in k much larger than crosses indicate) with a 
single Majorana fermion mode (arrow). Inset: zoom-out at higher 
energies. Conical structure at k = ±7r is due to the nodes in the 
day-wave order parameter. 

Dirac cone (red) close to k — ~tt (exact position is set by 
/ii) lives on the left hand side of the QSHI and is thus unaf- 
fected by the SC. The low-lying band at E > (green) lives 
on the unaffacted right hand side of the SC and thus gives 
the doubly degenerate zero energy states always present at the 
surface of dry-wave SCs. For t — 0, the zero energy states 
on the left hand side of the SC creates a mirror-symmetric 
band at E < to the green band. However, for finite t this 
band (blue, crosses) is found mainly at non-zero energies and 
it has only one single crossing of E = (arrow). This is 
a single Majorana fermion mode and its /c-value is the same 



as the fc-value for the (now gapped) Dirac apex of the QSHI 
edge state. However, this mode joins the nodal quasiparticle 
excitations arising from the bulk d-wave order parameter at 
k = ±7r. Thus, this single Majorana mode has no excitation 
gap to topologically trivial states. Thus, the single Majorana 
fermion at E = has no excitation gap to topologically triv- 
ial states. The band structure in Fig. |4]is qualitatively similar 
for all finite t and moderately small to large However, for 
Hi m the would-be Majorana mode merges with the conical 
d-wave nodal structure at k — ±ir and disappears, i.e. there 
is then no state at zero energy that is localized at the QSHI- 
SC interface. We also note that in order to achieve the band 
structure in Fig. [4] self-consistency is crucial. With a step- 
function order parameter Aj, the band for the zero energy 
states at the QSHI-SC interface is essentially mirror symmet- 
ric to the equivalent band residing on the unaffected side of the 
SC, and thus no single Majorana mode is created. This is in 
sharp contrast to our results for a s-wave SC and is due to the 
intricate interplay between the QSHI edge state and the zero 
energy states of the d^y-wave SC. To summarize this section, 
we have shown that a single Majorana fermion exists only for 
doped QSHI-d^y-wave SC structures and that there is no ex- 
citation gap for this mode. Thus, from an experimental point 
of view, we do not expect this Majorana fermion to be stable 
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for quantum computing operations [ 

The author thanks Jacob Linder and Eddy Ardonne for very 
valuable discussions. 
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